Spin waves and magnetic interactions in LiCu2 02. 
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The quasi-one-dimensional helimagnet LiCu2 02 was studied by single crystal inelastic neutron 
scattering. The dispersion relation of spin wave excitations was measured in the vicinity of the 
principal magnetic Bragg reflection. A spin wave theoretical analysis of the data yields an estimate 
of the relevant exchange constants and explains the mechanism of geometric frustration that leads 
to helimagnetism. It is found that the simple antiferromagnetic Ji — J2 model that was previously 
proposed is inadequate for LiCu2 02. The experimental findings are generally in a qualitative agree- 
ment with first principal calculations of [A. A. Gippius et ai, Phys. Rev. B 70, 020406 (2004)], 
though certain important discrepancies remain to be explained. 

PACS numbers: 



I. INTRODUCTION 



LiCu2 02 is a quasi-one-dimensional magnetic material 
that has recently attracted a great deal of attentioniSi^i^ 
Muon spin resonance^, neutron scattering experiments^ 
and NMR studies^ clearly demonstrated that this S = 
1/2 system has an incommensurate helimagnetic ground 
state. In our previous work^^ we postulated a simple 
model of magnetic interactions that seemed to account 
well for the observed properties. It was proposed that 
LiCu202 features weakly coupled zig-zag 5 = 1/2 chain 
with competing nearest-neighbor Ji and next-nearest- 
neighbor J2 antiferromagnetic interactions. Due to ge- 
ometric frustration this model exhibits an unusual com- 
petition between a quantum-disordered gappedS, ground 
state with commensurate spin correlationsjS and semi- 
classical helimagnetic incommensurate state. However, 
independent first-principle numerical studie o^d° sug- 
gested a totally different model for the geometric frustra- 
tion in LiCu202. According to that work, the nearest- 
neighbor Ji interactions are actually ferromagnetic, and 
it is the unusually large fourth-nearest neighbor antifer- 
romagnetic coupling constant J4 that causes geometric 
frustration and favors helimagnetism. A similar scenario 
was recently proposed for a very similar isostructural he- 
limagnet NaCu202«ii 



Bulk magnetic susceptibility data for LiCu202 seem 
to be compatible with both models^iiS and are thus in- 
conclusive. The only direct way to determine (effective) 
exchange constants is by means of single crystal inelastic 
neutron spectroscopy. This technique can probe the dis- 
persion relation of spin wave excitations in the system. In 
the present paper we describe the results of such a study. 
We find that the actual topology of magnetic coupling is 
indeed more complex than originally envisioned. 



Structural considerations 



LiCu2 02 has a layered charge-ordered orthorhombic 
crystal structure (space group Pnma, a = 5.730(1)1, b = 
2.8606(4)1^ and c = 12.417(2)1), as described in detail 
m Refs.a^ The magnetism is due to double chains of 
Cu^"*" ions that run along the crystallographic b axis. The 
period of each "leg" of the double spin chains is equal to 
b. The two legs are offset by b/2 relative to each other and 
thus form a "triangular ladder", as illustrated in Fig.Q^. 
These double chains are arranged in layers parallel to 
the (a, b) plane and are effectively separated along the c 
direction by planes of non-magnetic Cu'*' ions. 

The approximate magnetic structure was determined 
in Ref.lgand is illustrated in Fig. 1 of that paper. The or- 
dering temperature is Tc = 24 K. Each double spin chain 
carries a planar helix of spins. The magnetic propaga- 
tion vector is qo = (0.5, 0) with ( ss 0.827. Consecutive 
spins on each rung are almost parallel, being rotated rela- 
tive to each other by an angle (j) = 27rC. Within each dou- 
ble chain any nearest-neighbor spins from opposite legs 
are almost antiparallel and form an angle (j)/2 — tt(. Pre- 
liminary neutron diffraction experiments were consistent 
with (a, b) being the spin rotation plane, but independent 
NMR studies clearly show an out-of-plane component ii 
This suggests that the spin rotation planes are, in fact, 
tilted relative to a high-symmetry orientation. The issue 
deserves further investigation, perhaps involving more ac- 
curate diffraction measurements. However, if magnetic 
anisotropy effects are negligible compared to Heisenberg 
exchange interactions (a reasonable assumption for Cu^"*" 
ions), the actual arrangement of spin rotation planes 
should not affect spin wave dispersion relations that are 
the focus of the present study. 

The main goal of the inelastic neutron scattering ex- 
periments described below is to establish the hierarchy of 
exchange interactions. The corresponding coupling con- 
stants are schematically shown in Fig. and define a 
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FIG. 1; (a) A schematic view of exchange interactions be- 
tween magnetic Cu^^ ions in LiCu2 02. (b) An equivalent 
Bravais lattice of spins obtained by displacing every other 
Cu^^ ion in the original non-Bravais lattice by the vector d. 
It is assumed that the spin Hamiltonian remains intact upon 
this transformation. 

model Heisenberg Hamiltonian for LiCu202: 

id 

+ J4Si,j^i+4J + J±^i,jSi,j + l] , (1) 

where the index i labels consecutive spins in each double 
chain, as shown in Fig.^, and j labels the double chains. 

II. EXPERIMENTAL 

In the present study we employed a 5 g single-crystal 
sample of LiCu202. Crystal mosaic was not particularly 
good: irregular, with as much as 4° FWHM spread. This 
circumstance imposed certain constraints as will be dis- 
cussed below. The measurements were carried out in two 
separate series of experiments. Neutrons with a fixed fi- 
nal energy oi Ef — 13.5 meV were used at the HB-1 
thermal 3-axis spectrometer at the High Flux Isotope Re- 
actor at ORNL (Setup 1). Pyrolitic graphite PG (002) 



reflections were used for monochromator and analyzer. 
A PG filter was installed after the sample to eliminate 
higher-order beam contamination. The setup employed 
48' -80' -80' -240 collimators. Alternatively, we utilized 
cold neutrons with a final energy fixed &i Ef — b meV at 
the TASP 3-axis instrument installed at the SINQ spal- 
lation source, Villigen PSI, Switzerland (Setup 2). PG 
monochromator and analyzer were used with no addi- 
tional coUimation and a cold Be filter positioned after 
the sample. In all cases the sample was mounted with 
the (a, h) plane coinciding with the horizontal scattering 
plane of the instrument. Sample environment was a stan- 
dard "Orange" -type flow cryostat. Most of the data were 
taken at T = 1.5 K. 

Data collection procedures were aimed at minimizing 
the adverse effects of the broad mosaic spread of the sam- 
ple. A large mosaic translated to a poor wave vector 
resolution perpendicular to the momentum transfer q. 
Since the spin wave dispersion is steepest along 6*, the 
best focusing conditions are obtained for q pointing close 
to that direction. Considering the intensity reduction 
due to the effect of magnetic form factors at large 
an optimal "window" for low-energy spin wave measure- 
ments was identified in the vicinity of the (0.5,0.827,0) 
and (0.5,1.173,0) magnetic Bragg peaks. Another ad- 
vantage of this reciprocal-space region is that it is rel- 
atively clear of spurious scattering and phonons, that 
were carefully checked for. Most scans were repeated 
at r = 75 K, i.e., well above Tc, to verify that the sig- 
nal observed at low temperature is indeed of magnetic 
origin. In some cases smoothed high-temperature scans 
were used as background for point-by-point subtraction 
from the corresponding low-temperature data sets. 

Another important technical problem that had to be 
dealt with is crystal twining. In LiCu202 twinning oc- 
curs at a microscopic leveli^ The a axis of one type of do- 
mains coincides with the h axis of the other domain type, 
since a ~ 2b. The immediate consequence is that neu- 
tron diffraction and inelastic scattering necessarily detect 
a joint signal originating from both domain types. As will 
be discussed in detail below, the contributions from dif- 
ferent domains and different spin wave branches within 
each domain could be reliably separated only in the direct 
proximity of the (0.5,0.827,0) and (0.5,1.173,0) peaks. 
Even so, twinning reduces the effective sample volume by 
a factor of two. To avoid any confusion we will use the 
indexes h, k and I to label wave vector components in 
the coordinate system defined by the reciprocal lattice of 
one particular domain type (type- A). Indexes hs = 2fc, 
ks — h/2 and Is = I will refer to the other crystallo- 
graphic domain (type-B). 



III. RESULTS 

A typical constant-energy scan collected using Setup 
1 is shown in Fig. |21 The non-magnetic background for 
this scan (featureless and typically about 80 counts per 5 
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FIG. 2: A constant-_E scan along the (—0.5, k, 0) reciprocal- 
space rod measured in LiCu2 02 at T = 1.5 K at an energy 
transfer of 5 meV (symbols) using Setup 1. The lines are a 
simulation based on the measured spin wave dispersion rela- 
tion and the known resolution function of the instrument. 



min.) was measured at T = 75 K. A linear fit to the mea- 
sured background was subtracted from the data shown. 
The four prominent peaks seen in this scan can be at- 
tributed to an acoustic spin wave emanating from two 
magnetic Bragg peaks at (0.5, 0.827, 0) and (0.5, 1.173, 0). 
Representative constant-g scans collected using the same 
setup are shown in Fig. |31 where a flat background has 
been subtracted from the data. All inelastic peaks ob- 
served for 0.75 ^ k < 1.25 and Hcu < 8 meV can be 
associated with a single branch of excitations. The cor- 
responding dispersion relation along the crystallographic 
b axis was extracted from the measured scans using Gaus- 
sian fits. The result is plotted in circles in Fig.^, where 
solid and open symbols indicate constant-q and constant- 
E measurements, respectively. Any attempts to follow 
the observed spin wave branch to higher energies were 
not successful. The mode's intensity drops progressively, 
while the background increases and becomes structured. 
As will be explained below, part of this problem may 
be due to a multitude of additional spin wave branches 
from both types of domains. The situation was fur- 
ther aggravated by limited resolution, phonon scattering 
and instrument-related spurious peaks. For this reason, 
throughout this paper we shall limit the discussion to ex- 
perimental data collected at low energies in direct prox- 
imity of the two above-mentioned magnetic Bragg peaks. 

Additional measurements were performed using Setup 
2 to determine the dispersion along the a* axis. Typi- 
cal scans and the measured dispersion curve are shown 
in Figs. [5l and 2|3, respectively. The important zone- 
boundary scan at (0, 0.827, 0) was also measured using 
Setup 1 in the range 2-15 meV and found to be fully 
consistent with that shown in Fig. [S]:. 
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FIG. 3: Typical constant-q scans measured in LiCu202 (sym- 
bols) at T = 1.5 K. The lines are as in Fig. H The shaded 
area in the lower scan is an "accidental Bragg" spurious peak 
originating from 2ki scattering in the monochromator, (0, 2, 0) 
Bragg scattering in the sample and inelastic thermal-diffuse 
scattering in the analyzer. The spurious peak appears much 
narrower than the experimental resolution (about 2 meV 
FWHM). 



Data analysis 

In order to extract the relevant exchange constants 
from the measured dispersion curves the data were ana- 
lyzed in the framework of semiclassical spin wave theory 
(SWT)pia 



1. Equivalent Bravais lattice 

At a first glance, calculating the dispersion relation is 
quite tedious, since we are dealing with a non-Bravais lat- 
tice of spins with four magnetic ions per unit cell. Within 
the approximation of the Hamiltonian (^), however, any 
interactions between the double-chain layers are ignored, 
and only two spins per unit cell remain. Moreover, the 
spin network can be made equivalent to one on a Bra- 
vais lattice by displacing the atoms as shown in Fig. 
The dynamic structure factor S'o(q, of this "straight- 
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FIG. 4: Spin wave dispersion measured in LiCu2 02 at 
T = 1.5 K (symbols). Shaded and open circles are data points 
obtained from constant--B and constant-g scans, respectively. 
The lines are a fit to the data, as described in the text. Heavy 
solid, heavy dashed and thin dotted lines correspond to Mod- 
els 1, 2 and 3, respectively. In both panels the ellipses rep- 
resent the FWHM of the instrument resolution function in 
the appropriate projection (open ellipse) and section (shaded 
ellipse) . 

ened out" spin network is directly related to that of the 
original model through: 

S{ci,uj) = 5o(q,c^)cos2^ + 

+ 5o(q+(0,l,0),c^)sin2^, (2a) 



qd = 27r/i(5^ + 2t:15:,, (2b) 

where q = {h,k, I) is the wave vector transfer and d is the 
relative displacement of the two legs in each double chain: 
d = aSx + c6z- For the crystal structure of LiCu2 02 
we have 5^ = 0.24 and 5^ = 0.19. It is important to 
emphasize that Eq. [3 is exact: it does not imply any 
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FIG. 5: Typical scans measured in LiCu202 using Setup 2 
(symbols) . The lines are as in Fig. |21 

particular properties of S'^"^ (q, uj) and ^(q, lj), and rehes 
only on the definitions of these correlation functions. 

2. SWT dynamic structure factor 

A very useful recipe for calculating the SWT dynamic 
structure factor of a simple helimagnet on a Bravais lat- 
tice can be found in Ref. The magnon dispersion 
relation is given by: 

iol = Al-Bl (3a) 

= 25{ Y + \ [-/q-qo + 'Vqo] " ^qo } > (3b) 
Bq = 25{-^ - - [Jq_qo + Jq+qJ}, (3c) 

where, for our model Hamiltonian, 

Jh.k,i — Ji cos(7rfc) + J2 cos(27rfc) -|- 

+ J4 cos(47rfc) + Jj_ cos(27r/i). (4) 

The requirement that qo minimizes J(q) fixes the rela- 
tion between Ji, J2, J4 and C- 

Ji + 4 J2 cos(7rC) -I- 16 J4 cos(7rC) cos(27rC) = 0. (5) 
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TABLE I: Sets of exchange constants obtained by fitting tfie calculated spin wave dispersion relation to the experimental 
curves in comparison with first-principle calculations. 





Model 1 


Model 2 


Model 3 


LDA" 


Jl 


6.4(0.95) meV 


105.5 (8.0) meV 


(fixed) 


0.4 meV 


J 2' 


-11.9 meV 


33.8 meV 


-14.0 meV 


-8.1 meV 


Ji 


7.4 (0.5) meV 


-1.6 (0.2) meV 


7.5 (0.1) meV 


14.4 meV 


J± 


1.8(0.2) meV 


0.23 (0.02) meV 


6.8 (0.4) meV 


5.7 mcV 



"From Ref. 7 

''Not refined: fixed by Eq.lU 



Conveniently, J± does not enter this expression. The 
single-magnon cross section has three contributions: 



'S'o"(q,^) 



S{uJ-Uq), 



(6a) 



l ^q-qo B^-'^o siu^-u;^^^,), (6b) 



8 ^q+qo 



6{UJ ~UJq+q„). (6C) 



Here S'-'-(q, w) represents fluctuations of the spin com- 
ponent perpendicular to the plane of the hcUx, while 
S~^{q,uj) and S^{q,uj) are in-plane excitations. The to- 
tal SWT cross section for LiCu202 is related to 5o(q, ti-") 
through Eq. [3 

It now becomes clear why measuring spin waves in 
LiCu202 is so technically challenging: there are six spin 
wave branches in each crystallographic domain. As a re- 
sult, in any experiment one has to deal with a total of 
twelve spin wave branches, all of which are acoustic and 
therefore contribute to scattering at low energies. Also, 
compared to a simple antiferromagnet on a Bravais lat- 
tice, the scattering by each individual branch is typically 
much weaker, since intensity is effectively re-distributed 
between twelve modes. The existence of a suitable, albeit 
narrow, measurement window that was exploited in our 
experiments is a fortunate coincidence. 



Fits to experimental data 

The energy of the spin wave branch observed in our ex- 
periments clearly goes to zero at the principal magnetic 
Bragg peaks of the A-type domain. It must therefore be 
associated with the S-^ (q, oj) mode and the first term in 
Eq. I2al An analysis of the polarization of spin wave exci- 
tations is beyond the scope of this work. Nevertheless, we 
can point out that (a, b) being the spin rotation plane^ 
is consistent with a strong S'^(q, w). Indeed, the cor- 
responding polarization factor for unpolarized neutrons 
ia a maximum for in-plane momentum transfers^ As a 
first step in our analysis we ignored all other spin wave 
branches and used Ea. l3al to fit the experimental data. It 
is easy to show that there are exactly two sets of exchange 
parameters that exactly reproduce (i) the experimentally 



determined incommensurability parameter C, (ii) the ex- 
citation energies at the zone boundaries (0.5,1,0) and 
(0, 0.827, 0), and (iii) the spin wave velocity along the b* 
direction near the magnetic Bragg peaks. Using these 
two sets of parameters as initial points, we employed 
a least-squares algorithm to best-fit the data under the 
rigid constraint set by Eq. [S] The two resulting opti- 
mized sets of parameters are listed in Tableland will be 
referred to as Models 1 and 2, respectively. For a pur- 
pose that will be made clear in the next section, we also 
attempted a fit to the data taken only along the b-axis 
while fixing Ji = 0. The measured a-axis dispersion was 
not included in this fit. The resulting set of exchange 
constants is also listed in Table |2 and will be referred 
to as Model 3. Dispersion relations for the out-of-plane 
branch and the first term in Eg . I2al calculated for type- A 
domains using Models 1 through 3 are plotted in lines in 
Fig.H 

The next step is to understand why the other 11 spin 
wave branches were not observed in the studied wave vec- 
tor and energy ranges. All twelve dispersion curves calcu- 
lated using Model 1 are plotted in Fig.|S^. The intensities 
of each branch were scaled by the energy transfer (to get 
rid of the l/uj factors in Eq.|HJl and are plotted in Fig.lHJ). 
The neutron polarization factor was not included in this 
calculation, and the maximum value (unity) was assumed 
for each branch. One readily sees that all branches orig- 
inating from type-B domains and the 1st term in Ea. l2al 
are extremely weak in our area of interest. The same is 
true for all spin waves associated with type-A domains 
and the 2nd term in Eq. I2al Contributions from type- 
B domains and the 2nd term in Eq. I2al are either weak 
or at very high energy. In the end, only the "princi- 
pal" modes originating from type-B domains and the 1st 
term in Ea. l2al are relevant to our experiments. Of these 
three branches the two in-plane modes are almost en- 
tirely outside the studied energy range and do not affect 
any of the scans collected. As a result, only one out- 
of-plane mode is seen. Similar arguments can be made 
regarding the dispersion measured along the a-axis. The 
simulation also explains why we were unable to follow 
the "principal" out-of-plane mode to higher energies. At 
around 8 meV energy transfer any inelastic scans become 
crowded with a host of peaks associated with other spin 
wave branches that can not be reliably resolved due to 
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FIG. 6: (a) Dispersion of all twelve spin wave branches along 
the (0.5, fc, 0) reciprocal-space rod in twinned LiCu202 crys- 
tals calculated using Model 1. (b) The corresponding struc- 
ture factors scaled by energy transfer. In both panels the out- 
of-plane modes are plotted in red. Modes shown in blue and 
green lines are polarized in the plane of spin rotation. Thick 
and thin lines of all types correspond to the first and second 
terms in Eg. I2al respectively. Solid and dashed lines refer to 
type-A and type-B crystallographic domains. Symbols are as 
in the top panel of Fig. 2] 



limited experimental resolution. 

From Fig.0|it is apparent that the studied energy and 
wave vector transfer range the dispersion relations pro- 
duced by Model 2 are virtually identical to those result- 
ing from Model 1. Both models fit the available exper- 
imental data rather well. The only small discrepancy 
is seen along the a-direction where the observed spin 
wave velocity near q — (—0.5,0.827,0) appears slightly 
larger than calculated, while the zone-boundary energy 
at q = (0,0.827,0) is slightly smaller. This discrepancy 
is easily explained if one assumes a small anisotropy- 
induced gap A in the spin wave spectrum at the magnetic 
Bragg peak position. From the existing scans it is pos- 
sible to tentatively estimate A « 1.5 meV, though more 



high-resolution data will be needed for a reliable determi- 
nation. As far as only the &-axis dispersion is concerned. 
Model 3 seems to work just as well as the other two sets 
of parameters, and for the (0.5, k, 0) reciprocal-space rod 
produces almost identical curves. However, this good 
agreement under the constraint Ji — necessarily re- 
sults in a very poor agreement with the a-axis dispersion 
relation. The corresponding bandwidth calculated us- 
ing Model 3 (thin dotted line in Fig. ^p) is considerably 
larger than the observed value. In other words, assuming 
Ji = in Hamiltonianncan not be made consistent with 
the bulk of the experimental data. 

The resulting parameter values for Model 1 were used 
to simulate the measured inelastic neutron scans. For 
this purpose we utilized the magnetic dynamic structure 
factor as given by Eq. |23 This model cross section was 
scaled by the magnetic form factor for Cu^^ and nu- 
merically convoluted with the 4-dimensional spectrom- 
eter resolution function. The latter was calculated in 
Gaussian form using the Cooper-Nathans approximation. 
The evolution of the corresponding FWHM resolution el- 
lipsoid in the surveyed areas of {E, q) space is shown in 
Fig. 21 An overall scaling factor was chosen to best-fit 
the measured scans. The results of such simulations are 
shown in solid lines in Figs.l^HS Overall, the experimen- 
tal peak shapes are well explained by resolution effects 
alone. 



IV. DISCUSSION 

As mentioned in the introduction, the knowledge of 
spin wave dispersion in an ordered magnetic material 
can, in principle, provide comprehensive information on 
the exchange interactions. This requires measurements 
of dispersion curves across the entire Brillouin zone. In 
the particular case of LiCu202, twinning, less than per- 
fect sample mosaic, and a complex spectrum limit us to 
only observing spin waves in the direct proximity of the 
principal magnetic Bragg reflections. For this reason the 
resulting sets of parameters can be regarded as unique 
only under the assumption of the validity of Hamiltonian 
Hand the SWT approximation. 

The Hamiltonian has a solid justification in the crys- 
tal structure, and is well supported by first-principle cal- 
culations. However, a vital question is to what extent 
an analysis based on essentially classical spin wave the- 
ory can be valid in the case of a quasi-low-dimensional 
frustrated magnet such as LiCu2 02? At the very least, 
the exchange constants obtained from such an analysis 
will be renormalized compared to their actual values due 
to quantum corrections. Nevertheless, even for highly 
anisotropic or frustrated systems, acoustic modes can be 
associated with waves in the order parameter field. In 
this case a renormalized SWT with some effective set of 
exchange parameters typically works rather well. Just 
two examples are the weakly-coupled 5' = 1/2 chains 
compound KCuFa (Ref. 14) and the quasi-2D frustrated 
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helimagnet Cs2CuCl4."'^ - Deviations from classical behav- 
ior in LiCu2 02 do not appear particularly severe. The 
typical energy scale of exchange interactions computed 
by LDA exceeds the temperature of 3D long-range or- 
dering by less than an order of magnitude. While the 
absolute value of the ordered moment has not been de- 
termined to date, that in the very similar helimagnet 
NaCu202 is quite large (0.64 /ib )-^^ We therefore have 
good reason to assume that renormalized SWT will work 
well for our system. 

Turning now to the two solutions that we obtained we 
note that Model 2 almost exactly corresponds to the sim- 
ple Ji — J2 model that we originally proposedn^ However, 
unless the renormalization of exchange parameters is se- 
vere and all actual in-chain exchange constants are much 
smaller than the effective SWT ones, the estimated pa- 
rameters are not consistent with a previous analysis of 
magnetic susceptibility.^ We thus conclude that Model 1 
is a much more likely to be an adequate description of 
LiCu202. It's key features, namely a ferromagnetic J2 
bond and a substantial J4 coupling constant, are sim- 
ilar to those of the LDA calculations of Ref. 0. Geo- 
metric frustration in LiCu202 is thus similar to that in 
NaCu202ii 

The main difference between the exchange parameters 
in Model 1 and those emerging from LDA calculations 
of Ref. is in the magnitude of Ji. As illustrated by 
Model 3, assuming that Ji is negligibly small makes the 
SWT calculation incompatible with the measurements. 
We conclude that either LDA severely underestimate Ji , 
or this exchange constant is renormalized by over an or- 
der of magnitude in the SWT. Of course, one can never 
entirely rule out the possibility that both the LDA cal- 



culations and our SWT analysis are missing some cru- 
cial terms in the Hamiltonian. Further experimental and 
theoretical insight will be required to fully resolve this 
remaining mystery. 



V. CONCLUSION 

In summary, based on inelastic neutron scattering mea- 
surements we can conclude that the simple J1-J2 model 
that we originally proposed does not apply to LiCu202. 
The frustration mechanism is more complex and resem- 
bles that proposed in Refs. F^IToL It involves a competi- 
tion between a combination of antiferromagnetic Ji and 
ferromagnetic J2 interactions against an additional anti- 
ferromagnetic long-range J4 coupling. We find the three 
corresponding exchange constants to be of comparable 
absolute strength. This discrepancy between an SWT- 
based interpretation of the experimental data and first- 
principle calculations remains to be explained. 
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